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Asymmetry of Superconductivity in Hole- and Electron-doped Cuprates: 
Explanation within Two-Particle Self-Consistent Analysis for the Three-Band Model 
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In the hole-doped cuprate superconductors, the superconducting transition temperature T c ex¬ 
hibits a dome-like feature against the doping rate. By contrast, recent experiments reveal that T c in 
the electron-doped systems monotonically increases as the doping is reduced, at least up to a very 
small doping rate. Here we show that this asymmetry is reproduced by performing a two-particle 
self-consistent analysis for the three-band model of the C 11 O 2 plane. This is explained as a combined 
effect of the intrinsic electron-hole asymmetry in systems comprising Cu3d and 02 p orbitals and 
the band-filling-dependent vertex correction. 

PACS numbers: 74.25.Dw, 74.72.-h, 74.20.Pq 


Despite the long history, there still remain various 
unsolved problems in the study of the high-T c cuprate 
superconductors. The striking difference in the dop¬ 
ing dependence of the superconducting transition tem¬ 
perature T c between the hole- and the electron-doped 
materials is among those unresolved issues. It is well 
known that in the hole-doped case, T c exhibits a dome¬ 
like feature against the doping rate, namely, T c first in¬ 
creases upon doping (underdoped), then yields a maxi¬ 
mum value (optimal), and finally decreases with further 
doping (overdoped). On the other hand, it was known 
for the electron-doped cases that T c abruptly appears as 
soon as the antiferromagnetism is lost with doping, and 
monotonically decreases as the doping rate increases. Re¬ 
cent experiments show that the antiferromagnetism can 
be suppressed down to very small doping rate, or even 
in the mother compound, when the apical oxygens are 
ideally removed in the T'-type crystal structure of the 
electron doped cuprates. Then, it has been revealed that 
T c monotonically increases with decreasing the electron 
doping at least up to a very small doping rate (less than 
5 percent), and is suggested to be superconducting even 
in the non-doped mother compound 1 4 . 

There have been some theoretical studies of the dop¬ 
ing dependence of T c . The fluctuation exchange (FLEX) 
approximation 5 ! for the single band Hubbard model gives 
a monotonic doping dependence of T c ^, and therefore 
has difficulties in understanding the doping dependence 
of T c in the hole-doped cuprates. Some studies consid¬ 
ered superconducting fluctuation in FLEX to circumvent 
this problemP^. There have also been some studies that 
adopt methods capable of dealing with the strong corre¬ 
lation effect^® 5 . In some of those studies, T c exhibits 
a dome-like doping dependence, but in those cases there 
would be difficulties in understanding the recent experi¬ 
mental results for the electron-doped case. The electron- 
hole asymmetry of T c was studied in a two band model 
that explicitly considers the 02 p orbital, but there, the 
antiferromagetic phase was obtained in a wide electron 
doping range^ in contradiction to the experiments men¬ 
tioned above 2 ®. 

It has been suggested that the difference in the charac¬ 


ter of the mother compound (Mott insulator or not) be¬ 
tween the hole-doped and the electron-doped systems can 
be attributed to the difference in the electr onic structure 
originating from the crystal structure 4 12 17 18 . Namely, 
while the crystal structure of the single-layer hole-doped 
cuprates is composed of Cu -0 octahedra (T-type), that of 
the electron-doped cuprates is composed of Cu -0 squares 
(T'-type) and (ideally) has no apical oxygens. Due to 
this difference, the copper 3 d -oxygen 2 p level offset in 
the T'-type structure tends to be smaller than that in 
the T-structure. Since the d-p level offset is small in 
the electron-doped system,the on-site effective C 7 , when 
mapped to the single-band Hubbard model, is also small. 
One might expect that this difference in the crystal struc¬ 
ture, and hence the difference in the effective on-site U, 
can provide an explanation for the electron-hole doping 
asymmetry of T c . However, the inner layers of multi¬ 
layered hole-doped cuprates also do not have apical oxy¬ 
gens and therefore have the same lattice structure as that 
of the electron doped cuprates. Still, it is known that T c 
exhibits a dome-like doping dependence even within the 
inner layers 19 . Therefore, it seems difficult to attribute 
the electron-hole asymmetry of the doping dependence 
of T c to the absence/presence of the apical oxygens. The 
aim of the present study is to understand the origin of 
this electron-hole asymmetry. Here, we stress that in the 
present study we focus only on the (non-)monotonicity 
of the doping dependence of T c , and leave the issue of 
the metallicity or Mottness of the mother compound to 
future studies. 

We start by demonstrating that this electron-hole dop¬ 
ing asymmetry of T c is difficult to understand within 
the single band Hubbard model even when realistic band 
structures are considered. We perform first principles 
band calculation of HgBa2Cu04 (a hole-doped system) 
and Nd2Cu04 (an electron-doped system), and obtain 
tight-binding models constructing maximally-localized 
Wannier baskP 0 SI. Instead of the typical T-type hole 
doped system La2Cu04, we adopt HgBa2Cu04 because 
(i) it is known that the hybridization of the d z 2 orbital 
cannot be neglected in La2Cu04^, and (ii) the band 
structures of HgBa2Cu04 and Nd2Cu04 are very simi- 
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TABLE I. Nearest (£ 1 ), second (£ 2 ) and third (£3) neighbor 
hopping integrals for the single-band models. 



ti(eV) 

ta(eV) 

ts(eV) 

Nd 2 Cu0 4 

-0.457 

0.0866 

-0.0865 

HgBa2Cu0 4 

-0.464 

0.0907 

-0.0842 


TABLE II. Hopping integrals and copper 3 d -oxygen 2 p level 
offset in the three-band models. 


idp(eV) 

tpxPy ( e V) 

tpxPx ( e V) 

Adp(eV) 

Nd 2 Cu0 4 

1.18 

-0.621 

0.137 

1.83 

HgBa2Cu0 4 

1.26 

-0.632 

0.133 

2.06 

La2Cu0 4 

1.38 

-0.616 

0.0899 

2.73 


lar, so that we can concentrate purely on the electron- 
hole asymmetry. To take into account the electron cor¬ 
relation effect beyond those taken into account in the 
LDA/GGA level, the on-site interaction has to be treated 
by some many-body technique as has been done in pre¬ 
vious studies 25 27 . In the present study, we adopt the 
two-particle self-consistent method (TPSC) proposed by 
Vilk and Tremblay 2 ^. In this method, the interaction ver¬ 
tices in the charge and spin channel are approximated as 
different constants, and these constants are determined 
so that the correlation functions satisfy their sum rules 
that originates from the Pauli’s principle. It has been 
shown in refP 4 that TPSC gives a dome-like doping de¬ 
pendence of T c for the single band Hubbard model with 
nearest neighbor hopping only. 

The obtained hopping integrals for the single-band 
models are given in Table [TJ and the corresponding band 
structures are shown in Fig JIJ upper panels). The eigen¬ 
value A of the linearized Eliashberg equation for d -wave 
pairing, which is a measure of T c (see below), is shown 
against the band filling in Fig ( 2 ] For both the models 
of HgBa2Cu04 and Nd2Cu04, we set the on-site repul¬ 
sion as U/t = 8 and the temperature T/t = 0 . 08 , and we 
take 128 x 128 meshes and 4096 Matsubara frequencies. 
As shown in Fig| 2 j A varies monotonically in both the 
hole- and the electron-doped cases, namely, the dome¬ 
like T c variance against doping obtained for the nearest- 
neighbor-hopping-only case (inset of Figj 2 ]shows the dop¬ 
ing dependence of A for the £i-only model) is lost when 
a realistic band structure is adopted. 

Considering the previous studies mentioned in the in¬ 
troduction, it may be questionable whether we can repro¬ 
duce the experimentally observed electron-hole asymme¬ 
try within the single band model even if we take into 
account the electron correlation effects beyond TPSC. 
Namely, the doping dependence of T c would be either 
dome-like shaped or monotonic on both electron and 
hole-doped cases when the same values of U are taken. 
Hence, we now proceed to the three-band model that 
explicitly considers the in-plane oxygen 2 p XiV orbitals in 
addition to the copper 3d x 2_ y 2 orbitaP 9 . We first con¬ 
structed the five-band model composed of the copper 



FIG. 1. (color online). Band structure of the single-(upper 
panels) and three-band models(lower) for Nd 2 Cu 04 (left) and 
HgBa 2 Cu0 4 (right). 



band filling 


FIG. 2. (color online). Band filling dependence of the eigen¬ 
value A of the linearized Eliashberg equation for the single 
band models. For comparison, the inset shows the doping de¬ 
pendence of A for the single band model with nearest neighbor 
hopping only. 


3d x 2_ y 2 orbital and four in-plane oxygen 2 p XjV orbitals 
by using the maximally-localized Wannier basis. Sub¬ 
sequently, we obtained the three-band model by remov¬ 
ing two orbitals which are oriented in the direction 
perpendicular to the Cu -0 bond. The obtained model 
parameters and the band structure are shown in Table 
II] and Fig [I] respectively. For comparison, the model 
parameters for La2Cu04 are also shown. The parame¬ 
ter values of Hg and Nd systems can be considered as 
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quite similar, and especially the similarity of A d p can be 
noticed if we compare the values to that of the La sys¬ 
tem, which has smaller apical oxygen height compared 
to the Hg system. The similarity of Ad p between the 
electron-doped and the hole-doped materials is expected 
to become even more prominent if we consider multilayer 
hole-doped cuprates, where one or both of the apical oxy¬ 
gens are missing depending on the layer. This means 
similar values of the on-site U when mapped to single 
band Hubbard models, as mentioned in the introductory 
part. 

To analyze the three-band model, the TPSC approach 
should be generalized for multi-band system s. We follow 
the generalization of TPSC presented in Refs.® 51 . Let us 
briefly review TPSC for the multi-band Hubbard model. 
Hereafter we make use of the matrix form in the same 
way as Refs.™. 

The Hamiltonian of the three-band model is given as 


susceptibilities can be determined from the following sum 
rules derived from the Pauli principle: 

2 r p _ 

- Y = n P ~ 2 <*tyt n *4> . ( 5a ) 

q 

2J2 _ A 

Y =»v+ 2 K4 n /4> - ( 5b ) 

q 

where n y is the particle number per site of orbital 
/i, obtained from = n^. However the 

ansatz introduced here violates the electron-hole symme¬ 
try. Therefore if > 1, considering the electron-hole 
transformation, the ansatz should be modified as 


U sp 



(6) 


h = + A d P Y n d°( r ) 

r,r',cr r,cr 

+ Y'52 U v n A r ) n A r )’ (!) 

r n 

where c^ a (r) is a creation operator of an electron with 
spin a and orbital (i = d,p x ,p y at site r, n pa (r) = 
cj xcr (r)c AXCr (r) is a number operator, Ad p is the d-p level 
difference, U p is the on-site Coulomb interaction. The 
band filling n is defined as the average number of elec¬ 
trons per unit cell, so that n — 5 corresponds to the 
non-doped case. In the analysis for this model, we set 
the on-site interaction Ud = lOeV and U p = 5eV, tem¬ 
perature T = O.OleV. We employ 64x64 k -point meshes 
and 4096 Matsubara frequencies. 

In this three-band model, similar to the single-orbital 
case, the spin and charge susceptibilities are evaluated as 


X sp (q) = [1 - X°(q)U sp ] V(g), (2a) 

X° h (q ) = [1 + X°(q)U ch ] _1 X °(q), (2b) 

where x°(fc) is the irreducible susceptibility and L/ sp ( ch ) 
is the effective interaction matrix for the spin (charge) 
channel. The irreducible susceptibility is given by 

xSU(«) = -Jf E (k)G%(k + q), (3) 


using the bare Green’s function G°„(fc) = [(ie n + // — 
where p is the chemical potential and H(k ) 
is the matrix elements of the hopping term of the Hamil¬ 
tonian in the momentum representation. Here we abbre¬ 
viate the wave numbers and the Matsubara frequencies 
as k (for the fermionic case) or q (bosonic). 

Since we consider only U y as the interaction, introduc¬ 
ing the ansatz, 


U sp 


AMMLjT 

MM 


(4) 


where = 1 — n /J(T . Since n /; > 1 is satisfied for any 
band filling used in this study, we make use of the trans¬ 
formed ansatz. 

Using the obtained susceptibilities as described above, 
the dressed Green’s function G(k) is determined by 
Dyson equation: 

G(k)~ 1 = G^°\k)~ 1 — E(fc), (7) 

and the self-energy X(fc) is given by 

q 

+ UM h (q)U] lml , m ,G^ m ,(k- q ). (8) 

Solving linearized Eliashberg equation, 

AA W (k) = ]T MmAMGmM') 

k'rrii 

X A7722 7713 (k')G 77747773 ( k ), (9) 

the eigenvalue A and the anomalous self-energy A (k) are 
obtained. Here the singlet pairing interaction V s (q) is 
given by 

r s (q) = -U- | U s ^x sp (q)U + \u ch x ch {q)U , (10) 

where U is the interaction matrix for the bare vertex. 
The superconducting transition temperature T c is the 
temperature where A reaches unity. In the present study, 
we calculate A at a fixed temperature and use it as a 
measure for T c . 

Let us move on to the calculation results of the effective 
three-band model for HgE^CuCL and Nd 2 Cu 04 . The 
spin susceptibility = 0) an d the absolute 

value of the dressed Green’s function |Gdd(fe, ie n =o)\ are 
shown in Fig |3] and Figj4j respectively. As the number of 
electrons decrease from the electron-doped region (n > 5) 
to the hole-doped region (n < 5), peaks of the spin sus¬ 
ceptibility around ( 7 r, 7 r) and the T point are enhanced 
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FIG. 3. (color online). Spin susceptibility of the three-band 
models = 0).(a)Hg system, band filling n = 

4.85, (b)Hg system, band filling n = 5.0, (c)Nd system,band 
filling n = 5.0, (d)Nd system, band filling n — 5.15. 


FIG. 5. (color online). The doping dependence of the d -wave 
eigenvalue of the linearized Eliashberg equation A in the three- 
band model. The inset shows the doping dependence of A at 
a lower temperature T = 0.003eV. 



FIG. 4. (color online). Absolute value of the dressed Green’s 
function of the three-band models \Gdd(k,ie n =o)\. (a)Hg 
system, band filling n = 4.80, (b)Hg system, band filling 
n = 4.85, (c)Hg system, band filling n = 5.0, (d)Nd system, 
band filling n = 5.0, (e)Nd system, band filling n = 5.15. 


and the absolute value of the dressed Green’s function is 
suppressed. It can be seen in Figj3^a)(b) that the Green’s 
function is particularly suppressed around (7r, 0) (0,7r), 
namely the hot spots (see the figures in the supplemen¬ 
tary material, in which the hot spots are more clearly 
visible). 

These behaviors can be explained as a combined effect 
of the following two factors. First, since the d x 2_ y 2 or¬ 
bital is not half-filled due to the d-p hybridization, the 
d x 2_ y 2 orbital approaches the half-filling by decreasing 
the number of electrons. Because of this, to satisfy the 
sum rule for the spin susceptibility Xddddio) within the 
d orbital, the vertex U^ dd necessarily increases. There¬ 
fore the spin susceptibility increases with decreasing the 
number of electrons. Secondly, since the Fermi level 


approaches the van Hove singularity point of the band 
structure as the number of electrons is reduced, the spin 
susceptibility around the Y point is enhanced. The en¬ 
hancement of the spin fluctuation results in the increase 
of the self energy, which in turn suppresses the Green’s 
function. 

We show the band filling dependence of the d-wave 
pairing eigenvalue A of the linearized Eliashberg equa¬ 
tion (measure of T c ) in Fig (5] This result is consistent 
with the doping dependence of T c in both the electron- 
and the hole-doped region (except near the non-doped 
regime, which we will discuss later). As shown in the 
inset of Figj5j this feature remains at the temperature 
where the d-wave eigenvalue A is above unity near the op¬ 
timal doping rate (T = 0.003eV, 80x80 /c-point meshes, 
and 8192 Matsubara frequencies). This result can be 
interpreted as follows. Monotonic increase of T c in the 
electron-doped region as the number of electrons is re¬ 
duced arises from the enhancement of Xdddd^) ar °und 
(7r, 7r), which works in favor of the d-wave pair scatter¬ 
ing. As the band filling enters the hole-doped region, 
Xdddd^) further enhanced around ( 7 r, 7 r), so that A also 
increases. However, both the suppression of the Green’s 
function and the enhancement of Xdddd^) ar °und the Y 
point work against d-wave superconductivity, and there¬ 
fore A turns to decrease with further hole doping beyond 
5h = 0.15, where Sh = 5 — n is the hole doping rate. Thus, 
the doping dependence of the superconducting transition 
temperature is naturally understood in both the electron- 
and the hole-doped cases. 

Since the Mott transition is not described within the 
formalism used in the present study, in Fig(5]we show the 
result near n = 5 by dashed lines (the calculations have 
been done also at n = 5 nonetheless). The absence of 
the insulating state in the non-doped case is attributed 
to the insufficiency of the evaluation of the local electron 
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correlation effects. The inclusion of further electron cor¬ 
relation effects is left for future study. Nonetheless, we 
can expect that the inclusion of such effects will probably 
make the dome-like feature in the hole-doped region more 
prominent, while it should somewhat reduce the enhance¬ 
ment of A in the underdoped regime of the electron-doped 
case, which seems rather strong in the present calculation 
compared to experimental observations 23 . Hence, the 
inclusion of further local correction is likely to make the 
doping dependence of T c even closer to those observed 
experimentally. Another related issue is the pseudogap 
problem in the underdoped regime. This has been ad¬ 
dressed by TPSC in ref. 28 for the single band model, but 
the situation can be different in the case of the three-band 
model with realistic band structure. This also serves as 
an interesting future problem. 

To summarize, we have studied the doping depen¬ 


dence of superconductivity for the three band model 
of Nd 2 Cu 04 and HgBa 2 Cu 04 using the TPSC method. 
The eigenvalue of the Eliashberg equation A exhibits an 
optimal doping around n = 4.85 (hole concentration 
Sh = 0.15) in the hole-doped region and varies mono- 
tonically in the electron-doped region, consistent with 
the experiment. It is found to be understood naturally 
in terms of the electron-hole asymmetry due to the d-p 
hybridization and the band-filling-dependent vertex cor¬ 
rection. 

Part of the numerical calculations were performed at 
the facilities of the Supercomputer Center, Institute for 
Solid State Physics, University of Tokyo.This study has 
been supported by Grants-in-Aid for Scientific Research 
No.26247057 from the Japan Society for the Promotion 
of Science. 
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